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We present a general formulation of the time-dependent initial value problem for a quantum 
scalar field of arbitrary mass and curvature coupling in a Friedmann-Robertson- Walker (FRW) 
cosmological model. We introduce an adiabatic number basis which has the virtue that the 
divergent parts of the quantum expectation value of the energy-momentum tensor {Tat) are 
. . . isolated in the vacuum piece of (Tab), and may be removed using adiabatic subtraction. The 

O^ . resulting renormalized {Tab) is conserved, independent of the cutofi', and has a physically trans- 

parent, quasiclassical form in terms of the average number of created adiabatic 'particles'. By 
analyzing the evolution of the adiabatic particle number in de Sitter spacetime we exhibit the 
time structure of the particle creation process, which can be understood in terms of the time 
at which different momentum scales enter the horizon. A numerical scheme to compute (Tab) 
as a function of time with arbitrary adiabatic initial states (not necessarily de Sitter invariant) 
^•^ , is described. For minimally coupled, massless fields, at late times the renormalized (Tab) goes 

fT^ ' asymptotically to the de Sitter invariant state previously found by Allen and Folacci, and not 

to the zero mass limit of the Bunch-Davies vacuum. If the mass m and the curvature coupling ^ 
differ from zero, but satisfy m^ +^R — 0, the energy density and pressure of the scalar field grow 
linearly in cosmic time demonstrating that, at least in this case, backreaction efi'ects become 
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^vj ' significant and cannot be neglected in de Sitter spacetime. 
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I. INTRODUCTION 

Quantum particle creation effects are important in black hole physics [|l|^ , the generation of density pertur- 
bations and primordial gravitational waves , and reheating in inflationary cosmological models [Q , the damping 
of anisotropy in the early Universe ||5[] , and perhaps for even the very existence of a cosmological constant y,|7|] . 
In addition, the issue of quantum backreaction due to these effects could also be important for understanding 
black holes [|| and the dynamics of the early Universe g . 

In a curved spacetime background, it is not possible in general to define a unique vacuum state pffl. However, if 
the background expansion is sufficiently slow the notion of an adiabatic vacuum state defined via WKB expansions 
comes closest to the Minkowski vacuum |ll|. For FRW cosmological models, adiabatic regularization |13| has 
been shown to be particularly convenient, and equivalent to covariant point-splitting |13-n3|. 

In this paper we present a consistent, physically simple formulation of the general time-dependent initial value 
problem for a quantum scalar field of arbitrary mass m, and curvature coupling ^ in a background homogeneous, 
isotropic FRW cosmology, which is suitable for implementation on a computer. We present numerical results for 
particle creation in de Sitter space, and outline how the formalism can be applied to cosmological backreaction. 

We make use of an adiabatic particle number basis obtained via a time-dependent Bogoliubov transformation, 
in which the energy-momentum tensor of a quantum field in a FRW universe breaks up naturally into a 'particle' 
and a 'vacuum' contribution, the latter isolating the ultraviolet divergences. These divergences are then removed 
using adiabatic subtractions. The remaining finite energy density and pressure of the quantum field admit an 
intuitive quasiclassical interpretation in terms of the average number of particles present at a given time. These 
terms can be interpreted as the energy density and pressure of the produced 'particles' if the Universe entered an 
adiabatic era. This definition is motivated by and is a direct outgrowth of parallel definitions of particle number 
in time dependent problems in field theories in Minkowski space, such as scalar QED and A$'' theory [p6|-p^. 

It is well known that the concept of particle is ambiguous in time dependent backgrounds. Covariantly con- 
served quantities such as the renormalized expectation value of {Tab) can be defined without ambiguity (save the 
usual freedom associated with the additive finite renormalization of conserved local tensors), although the division 
of {Tab) into 'vacuum' and 'particle' components is not unique. However, if one requires that the intrinsically 
quantum ultraviolet power divergences of {Tab) arc isolated in the 'vacuum' component of {Tab), then the ambi- 
guity in the remainder is considerably reduced. In fact, this requirement is equivalent to matching the definition 
of 'vacuum' to an adiabatic basis up to at least second adiabatic order. One does not wish to match the vacuum 
to higher than fourth adiabatic order since this would introduce higher than fourth order time derivatives in the 
backreaction equations. Hence the only remaining ambiguity in the definition of the 'particles' resides in the 
treatment of the fourth adiabatic order terms, and is actually relatively mild. We introduce a precise definition 
of adiabatic basis in this paper which resolves the remaining ambiguity and use it to renormalize the energy- 
momentum tensor and study the particle creation process in de Sitter background. The resulting renormalized 
{Tab) has a physically transparent, quasiclassical form in terms of the average number of created 'particles', and 
is the appropriate source term for dynamical backreaction of the quantum field on the metric in semiclassical 
cosmology. Since our main interest here is the initial value problem in semiclassical gravity, we consider general 
FRW-invariant initial states. To prevent the occurrence of initial time singularities, we restrict attention to initial 
conditions which are also matched to the appropriate instantaneous fourth order adiabatic vacuum at the initial 
time |l|. 

The renormalized {Tab) we obtain is cutoff independent, and covariantly conserved. As a consequence of our 
analysis we demonstrate that these features do not depend on the choice of time parametrization {e.g., comoving 
versus conformal time), thus clarifying some previous discussions in the literature ^^. We have implemented a 
numerical scheme to compute the time evolution of the renormalized {Tab) which preserves covariant conservation 
and cutoff independence. Though our scheme is applicable to general FRW spacetimes, in this paper we have 
restricted ourselves to de Sitter space with closed spatial sections. We have checked covariant conservation and 
cutoff independence of the renormalized {Tab) for a variety of masses, couplings to curvature, spatial cutoffs, and 
non de Sitter- invariant initial states. 

We consider in some detail the time structure of the particle creation process for massive scalar particles in 
de Sitter space starting from a well-defined adiabatic initial state, which is not de Sitter invariant. The particle 
number in a particular momentum mode increases sharply when the physical wavelength of that mode crosses 
the de Sitter horizon. We show how the adiabatic particle basis is useful in interpreting the numerical results for 
the time evolution of the renormalized energy and pressure. We apply the same method to the massless case with 
^ = and ^ = 1/6. For the case ^ = 1/6, as expected from conformal invariance, there is no particle production, 
and the renormalized value of {Tab) is determined by the trace anomaly [Q. 



The case of M^ — m? + ^R = is important because the two hchcity states of gravitons obey a wave equation 
which corresponds to this value of the mass parameter M^, and there are several indications of strong infrared 
effects for such fields in de Sitter space. For example it is known that no de Sitter invariant state of the kind 
usually assumed, implicitly or explicitly, in inflationary models of the early universe exists for massless, minimally 
coupled scalar fields |p2| , and that the graviton propagator has unusual infrared behavior in de Sitter space p3| . 

We will see that the exact minimally coupled, massless case m? = ^ = is special in that in this case, at late 
times the renormalized {Tab) goes asymptotically to the constant de Sitter invariant value found by Allen and 
Folacci |22J , and not to the zero mass limit of the Bunch-Davies vacuum ||2l[ . For any other nonzero values of 
rn? and ^ such that Af ^ — the energy-momentum tensor grows without bound, linearly with cosmic time i, for 
adiabatic finite initial data. (We point out that this result is a quantum effect, and not a classical instability.) 
Hence at least in these cases it is clear that the backreaction of {Tab) on the metric cannot be neglected at late 
times. We have also investigated the case of a small mass, minimally coupled field, results for which will be 
reported elsewhere |24[] . 

The outline of the paper is as follows. In the next section we review scalar quantum field theory in general 
FRW spacetimes, which will serve to fix our notation and conventions. In Section III we introduce the adiabatic 
number basis by a certain time dependent Bogoliubov transformation on the exact mode functions. In Section 
IV we derive the energy density and pressure for a scalar field of arbitrary mass and curvature coupling in 
homogeneous, isotropic cosmologies as sums (or integrals) over the time dependent Fourier mode amplitudes and 
express the energy and pressure in the adiabatic number basis, showing explicitly how to isolate and remove 
the ultraviolet divergences in the vacuum contributions to these mode sums. In Section V we present numerical 
results for de Sitter space, for the massive, conformally coupled and minimally coupled scalars, for a variety of 
finite non de Sitter invariant initial states. In Section VI we present our conclusions and discuss the extension of 
our method to the full backreaction problem. A proof of the independence of the adiabatic method under time 
reparametrizations is given in the Appendix. 

II. SCALAR FIELD IN FRW SPACETIMES 

We consider spatially homogeneous and isotropic spacetimes which are described by the line element 

ds^ = -dt^ + g.jdx'dx^ = -dt^ + a^{t) dY? , (2.1) 

in the cosmic time i, comoving with freely falling observers. Here dS'^ is the line element of spatial sections with 
constant curvature, which can be labeled by k = —1,0, +1, viz., 
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dY?= '^'^ , + r^ (d^^ + sin^ g d(^' ) . (2.2) 

1 — Kr^ 

The conformal time variable ry is defined by the differential relation 

dt , , 

dr? = ^-r , 2.3 

a{t) 

and the function 

C{r,) ^ a" {tirj)) , (2.4) 



so that the line element (2.1) may be expressed in the equivalent form, 

ds^^C{r]) (-dry^+dE^). (2.5) 

It is the latter form that was considered by Bunch [E5[ , a reference to which some of our results in this and the 
next section may be compared. 

The time dependent Hubble function in comoving time is 



where the dot denotes differentiation with respect to comoving time t. The non-vanishing components of the 
Riemann, Ricci, and Einstein tensors, Gab ~ Rab — h^-Qab, are given by 
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Rtt ^ -'S{H + H^) , 

Rro = (H + -iH'^+'^] g,j 



Gtt = i H 



G,j=-l2H + 3H' + —]g,, 



Hence the the scalar curvature is 
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R^9"■"Rab = -Ru+9''R^,=& 



H{t) + 2H^{t) + ^ 



We will require the non-zero components of the conserved geometric tensor, 



Wi/.,^-L-^/v^d^.i?^ 

V^ Sg'"' J 



as well, which are given by 



2VaVbR - 2gabUR - -gabR^ + 2RRab , 



^^^Htt = -6HR + ii?2 _ g(-^2 ^ ^^^ ^ 

-2aR - 2HR + -i?2 - 2{H^ + H)R 
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(i)i/° = -eni? = &{R + 3HR) . 

In a conformally flat FRW spacetime the Weyl tensor, defined as follows 

Cabcd = Rabcd — 77 {QacRbd + QbdRac — OadRbc — 9bcRad) + -^ R (,9ac9bd — 9ad9bc) 

2 o 

vanishes. From this it follows that the geometric tensor 

2 1 1 

Hab = RaRcb — TyRRab — -zRcdR*^ 9ab + -J R 9ab , 



which satisfies 



V^ ^^'>Hab = 2R'-'WbClaa = 



is covariantly conserved in conformally flat FRW spacetimes as well. 

We consider in this paper a free scalar field <& with arbitrary mass and curvature coupling 
by the quadratic action. 



(2.7) 
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(2.9) 

(2.10) 

(2.11) 
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S=-\l V^ d^x [(Va$)g"^(V6$) + to2$2 ^ ^^^2] ^ 



in a general curved spacetime, where Va denotes the covariant derivative and g = det(ga6) 
for $ obtained by varying this action is 



(2.17) 
;, which is described 

(2.18) 
The wave equation 
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-ei?]$(t,x) = [-D + M^] $(f,x) ==0. 



(2.19) 



In a FRW spacetime expressed in comoving time coordinates (2.1) the D'Alembert wave operator on scalars 
becomes 



n$ = -5°''VaVf,$ = $ + 3Hi - 4jA3$ 



(2.20) 



where A3 is the Laplace-Beltrami operator in the three dimensional spacehke manifold of constant curvature k. 
This operator is diagonalized by the spatial harmonic functions Y]s_{x.) which satisfy 



- A3rk(x) = (fc2 - «;)rk(x) 



(2.21) 



where k = |k| takes on the discrete values of the positive integers 1,2, .. . in the case of k = 1, but becomes a 
continuous index over the non-negative real numbers in the case of k = 0, — 1. In the case of flat spatial sections, 
K = 0, the harmonic functions lk(x) become the ordinary Fourier plane wave modes e*'''^. In the compact 
case the harmonic functions are the spherical harmonics of the sphere S^, which are labeled by three integers 
k -^ (fc, £, m), the latter two of which refer to the familiar spherical harmonics on S'^ with £ < k—1. Thus a given 
eigenvalue of A3 labeled by k is fc^-fold degenerate. We note that in the present conventions the lowest mode is 
the k = 1 mode which is the constant mode on S^. The scalar spherical harmonics Ykim niay be normalized in 
the standard manner to satisfy 



S3 

when integrated over the unit three sphere, and 

fc-l £ 



d^^ ^fc%'rn'(x) Yklmi^) = Sk' kSf iS,n'm , 
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(2.22) 
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which is independent of x. 

Because of homogeneity and isotropy the wav e equation ( 2.19 ) separates in either comoving or conformal 
time coordinates, and the general solution of (|2.19) may be written in the form 



<i>(i,x) - a-i{t)J2 [ak/fc(t)>k(x) +aU-W>i:(x)] 



(2.24) 



where the time dependent mode functions fk{t) satisfy the ordinary differential equation, 

Mt) + nUt)Mt) = o, 

which is the equation of a harmonic oscillator with the time varying frequency, 
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(2.25) 



(2.26) 



We denote by w^ the first two terms in parentheses, which involve no time derivatives of the metri c, while Suj'^ 
denotes the remaining terms which are second order in time derivatives of the scale factor. In (2.24) we use the 
discrete notation or the sum ^j^, which is replaced by an integral over d^k in the cases k — 0,-1. 

Quantization of the scalar field is effected by requiring the equal time canonical commutation relation 



[$(i,x'),n(i,x)] =iM^(x',x) 



(2.27) 



where 11 = ^/—g^ = a'^$ is the momentum conjugate to $. Because of the completeness and orthonormality of 
the spatial harmonic functions ik(x) this canonical commutation relation is satisfied provided the creation and 
destruction operators obey 



[aicflk/] ^ "^kk' — Sk'kSpjpSm'r. 



(2.28) 



and the complex mode functions obey the Wronskian condition. 



hfl-hrk=-^^, (2.29) 



which they may be chosen to satisfy at the initial time to- From the mode equation of motion ( p. 25 ) this Wronskian 
condition will be fulfilled then for all subsequent times. 



We will consider states of the scalar field, which like the geometry (2.1), are also spatially homogeneous and 



isotropic. This implies that the expectation value of the particle number operator 

(a[ak> = (flkfll) -l = Nk, (2.30) 

can be a function only of the magnitude |k| = A;. In the context of Minkowski spacetime quantum field theory it 
has been shown that it is always possible to fix the bilinears [l^ ] 

{aUi) = (okflk) = , (2.31) 

by making use of the freedom in the initial phases of the mode functions fk , with no loss of generality. Precisely 
the same conclusion applies here. Because of the ubiquitous appearance of the Bose-Einstein factor in succeeding 
sections we will often make use of the the definition 

ak = 2Nk + l, (2.32) 

in what follows. The 'vacuum' state in this basis corresponds to the choice (7^ = 1. 

III. ADIABATIC NUMBER BASIS 



The observation underlying the introduction of the adiabatic basis is that the mode equation ( 2.25|) gener- 



ally possesses time dependent solutions which have no clear a priori physical meaning in terms of particles or 
antiparticles. The familiar notion that positive energy solutions of the wave equation correspond to particles 
while negative energy solutions correspond to antiparticles is quite meaningless in time dependent backgrounds, 
where the energy of individual particle/antiparticle modes is not conserved, and no such neat invariant separation 
into positive and negative energy solutions of the wave equation is possible. This is just a reflection of the fact 
that particle number does not correspond to a sharp operator which commutes with the Hamiltonian, i.e., parti- 
cle/antiparticle pairs are created or destroyed, and physical particle number is not conserved in time dependent 
backgrounds Pq] . 

It is equally clear physically that in the limit of slowly varying time dependent backgrounds there should 
be an appropriate slowly varying particle number operator which becomes exactly conserved in the limit of a{t) 
vanishing. Clearly th is slo wly varying particle number is not the Nk defined with respect to the time independent 
Heisenberg basis by ( p. 30 ) above. This N^ is part of the initial data, a strict constant of motion, no matter how 



rapidly varying the spacetime geometry is. The physical particle number at time t must be defined instead with 
respect to a time dependent basis which permits a semiclassical correspondence limit to ordinary positive energy 
plane wave solutions in the limit of slowly varying rifc(i). This adiabatic basis is specified by introducing the 
adiabatic mode functions fb6|j27f| 



A(0 = W^p^expHefe(i)) , with efe(t)=/ dt'Wkif), (3.1) 



which are guaranteed to satisfy the Wronskian condition ( 2.29 ) for any real, positive function Wk{t). The phase 



Qk{t) in ( p.l| ) can be measured from any convenient point such as the initial time io- 



If one requires fk{t) to be a solution of the exact mode equation ( 2.25 ), a second order non-linear equation 



for the frequency W'^ = il^ + j^ — j^ is obtained. However, the utility of the WKB form of this equation is 

that instead of solving for Wk exactly, it can be used to generate an asymptotic series in order of time derivatives 
of the background metric. Terminating this asymptotic series at a given order defines an adiabatic mode function 



fk{t) (to that order), which will no longer be exact solution of the original mode equation (2.25), but which 
can serve as a template adiabatic basis against which the exact mode functions fk{t) can be compared. This 
is accomplished by expressing the exact mode functions in the adiabatic basis by means of a time dependent 
Bogoliubov transformation of the general form P§,p3l ' 



fk{t) - ak{t)h{t) + Pk{t)fm ■ (3.2) 

The precise form of the time dependent BogoHubov coefficients ak and (3k will be determined once we specify 
completely the adiabatic basis functions fk, and these will be used to define a time dependent adiabatic particle 
number. 

The first few orders of the asymptotic adiabatic expansion of the frequency Wk in powers of the time variation 
of the metric are given by 
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wr--i^K+m' 



„,(2)2 __^2 , 3Wf^ T^r _o2 , 3u^' '^^■ 

Wt''=nl + ^^^^, etc. (3.3) 



Explicitly to second adiabatic order in the asymptotic expansion, 



where we have dropped third and higher order terms in the ellipsis. From this we evaluate 

-^"^Hri-^V%#(..2.H) (3.5) 



W^f V -^IJ 2C.2 

^k ^^k ^^k ^^k 

correct to adiabatic order three. Finally we have 

W^"^^ = WJf^ - ^^ ~ f ^ i?2 „ 1^—11 {R + 5RH + 2RH + 6RH'^) (3.6) 
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{5RH + 3RH + 19RH^) + -— ^ (h + 15HH + lOH^ + mHH^ + QQH^) 
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!^^ y- RH^ = (28HH + 19H^ + iMHH^ + b07H*) 

16 w^ 32 o;^ 

32 w^ 128cjii 

where terms of adiabatic order higher than four have been discarded. Notice that the terms fall off faster with large 
k as we develop the expansion to higher adiabatic orders, and therefore the adiabatic expansion is an asymptotic 
expansion in powers of 1/fc which matches the behavior of the exact mode functions in the far ultraviolet region. 
This is what makes the expansion useful for defining a particle number basis, since it can be used to isolate the 
divergences of the vacuum energy-momentum tensor. Since Tab is a dimension four operator in four spacetime 
dimensions, its divergences are obtained by expanding to adiabatic order four, and no higher orders are required. 



.(4) 



In fact all terms in W^ beyond the first three of (3.6) fall off faster than Wj, as fc ^ cx), and therefore will lead 



to convergent mode sums (integrals) and finite contributions to (Tab)- 

Our method will make use of these well-known results and take them one step further. Instead of demanding 
that fk(t) satisfy the exact mode equation ( ^.25 ) we will define Wk{t) to match the terms in W^ that fall off 



most slowly in k, in order to isolate the divergences of the energy-momentum tensor in the vacuum sector in 
the most convenient way. The exact mode functions fkit) will be written then as a linear combination of //c(i) 
defined by (^) and its complex conjugate. The fkit) defined in this way specify an adiabatic number basis in 
which particle creation in the homogeneous, isotropic state of the quantum field $ can be discussed in a physically 
meaningful way, free of ultraviolet divergences. 



Since the Bogoliubov transformation (3.2) may be viewed as a pair of canonical transformations in the phase 
space spanned by the field variables and their conjugate momenta, a complete specification of the transformation 
between bases requires a condition on the first time deriv ative s of the mode functions as well. The general form 
of this condition which preserves the Wronskian relation ( ^.29 ) is 



fk{t) 



-tWkit) 



Vk{t) 



akit)fk{t) 



tWk{t) 



Vk{t) 



(ik{t)im. 



(3.7) 



which introduces a second real function of time Vk ■ Usually this second independent function Vk , which contains 
only odd adiabatic orders, has been set equal to —Wk/Wk in earlier discussions of the adiabatic expansion. 
Hence the general adiabatic basis as a canonical transformation in phase space, where Wk and Vk are treated 
independently and on an equal footing, has not been fully realized. We will make use of the freedom to choose 
both Wk and Vk independently in our development. The precise choice of Wk and Vk will be postponed until 
the next section, where we present the detailed analysis of the ultraviolet divergences of the energy-momentum 
tensor. 



For any real Wk and Vk the hnear relations (3^) and (3/7) may be inverted to solve for the Bogoliubov 
parameters 



ak{t) = irkit) 

Pk{t) = -t.fk{t) 



.fkit)-[iWk{t) + ^]fk{t) 



fk{t) + tWk{t) - 



2 

Vk{t) 



hit) 



It is readily verified that these Bogoliubov parameters satisfy the relation 



(3.8) 



(3.9) 



for each k and for any choice of real Wk and Vk^ which is what is required for the transformation to be canonical. 
Making use of the mode equation (2.25) and the definition of the adiabatic functions in ( |3.l| ) we obtain the 
dynamical equations for the Bogoliubov coefficients 
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(3.10) 



An equivalent form of the canonical transformation in the Fock space of creation and destruction operators is 
given by |27| 



flk 



al(tya^(t) - P*k{t)al{t) , and al = ak{t)al{t) - l3k{tYa^{t) , 



(3.11) 



where the coefficients depend only on the absolute magnitude k = |k| by homogeneity and isotropy. We may now 
define the adiabatic particle number to be 



Nk{t) = (4(0ak(i)) = |afc|'(aj,ak) + |/3fc|'(akal) 

= (1 + |/3fe|2) Nk + |/3fe|2 (1 + Nk) = Nk + {l + 2Nk) |/3fe(t)p . 

If the initial condition on the mode functions is chosen to match the adiabatic functions, viz.. 



(3.12) 



fk{to) = fk{to) , and fk{to) 



'^Wk{t) + ^]fk{to)^ 



so that 



(3.13) 



Q!fc(io) = 1 I E^nd Pk{ta) = (adiabatic vacuum) 



then the last expression for the adiabatic particle number at arbitrary time t in ( 3.12 ) may be viewed as the sum 
of the particle number present at the initial time, Nk, plus the particles created by the time varying geometry. 



(3.14) 



|/3fcp, multiplied by the Bose enhancement factor, 1 + 2Nk = <Jk, to account for both spontaneous and induced 
particle creation processes. 



Because the rapidly varying phase Qk drops out of the definition of the particle number in (3.12), J\fk is a 
relatively slowly varying function of time. In fact, if Wk and Vk are chosen appropriately to match the adiabatic 
mode functions to a given order, Afk will be an adiabatic invariant to that order. On the other hand the bilinear 



Ckit) = (ak(t)ak(t)> = (1 + 2Nk)akP*k 



(3.15) 



is a very rapidly varying function of time, since the phase variable &k does not drop out of this combination. We 
can remove the leading part of this rapid phase variation by defining the quantities 



7^fc(t) = Re(Cfee-2^«'^) ^ a^ Re(afc/3^e-2«'=) , 
Jfc(t) EE Im(Cfee-2«'=) = afeIm(afe/3^e-2^®'=) , 



(3.16) 



which together with Mk span the space of three non-trivial real bilincars in the adiabatic creation and destruction 
operators. 

Because of the relation (p^) these three time dependent functions arc not independent, but rather are related 

by 
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(3.17) 



where ak is a str ict c onstant of motion, de pende nt only on the initial data. Differentiating (3.12) with respect to 
time and using (3.10), and the definitions (3.16), we obtain as well, 



—Uk ^nk( — + Vk 
dt \Wk 




^l - wl 



Vk , v^ 



(3.18) 



The two relations ( 3.17 ) and ( 3.18 ) imply that TZk and 2k may be eliminated in favor oi Afk and Afk, if desired. 
Explicitly, 



7^fe = \Ck\ cos(i?fe -Qk) ak^\Ck\ sin(t9fe - Cfe) 
with the angles dk and Qk defined by 

dk = sin"^ 



(3.19) 



Cfe = sin 
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(3.20) 



and \Ck\ is given in terms oi Afk and the constant Nk by ( |3.17 ) above. These relations are the general case of 
the S'C/(1, 1) structure of the group of Bogoliubov transformations characteristic of a Gaussian statistical density 
matrix, studied previously in the context of the leading order large N limit of flat space field theory |p6 17|. 

They will imply in the next section that all components of the energy-momentum tensor of a free scalar field in 
a general FRW background can be expressed in terms of the particle number Afk and its first time derivative, once 
the adiabatic particle basis is completely specified by the functions Wk and Vk- (This is a useful fact which suggests 
an extension of the semiclassical Boltzmann-Vlasov description of particle creation in cosmological spacetimes, 
along the lines of earlier results for particle creation in time dependent electromagnetic potentials p8|.) 

For completeness we give here also the first order differential equations obeyed by the functions TZk{t) and 
Tk{t), 
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In all of this discussion we are free to make any choice of the time dependent functions Wk and Vk in (3J_), (3.2), 
and (3.7) which is convenient for our purpose, since no physical quantity can depend on our choice of basis for 
the mode functions. The key point will be to use this freedom to choose the functions Wk and Vk so that the 
ultraviolet divergences in the energy density and pressure may be isolated in the vacuum sector and removed in a 
simple way. This will require that Wk and Vk be matched to the appropriate order of the corresponding adiabatic 
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functions W^ and —W^ l^k given for n = 2, and 4 by ( |3.4[) , and ( p.6[) , respectively. 



IV. THE ENERGY-MOMENTUM TENSOR 



The classical energy-momentum tenso r of a free scalar field in an arbitrary background gravitational field gab 
follows from variation of the action ( 2.18|) |ll|,p5[ , 



Tab = (Va$)(V6$) - ^g-'i(Vc$)(Vd$) - 2^Va($Vfa$) + 2^3,65='* Vd($V,$) 



and its trace is given by 



-5afc$ 



T = g'^^Tah = (6C - l)5°^Va($Vh$) - m^^^ _ $(_n + ^i? + m^)^ , 



where the last term vanishes by the equation of motion ( 2.19| ) . 
For a perfect fluid in a FRW spacetime we have 

Tab = V9ab + (e + p)WaWfc ■ 



(4.1) 



(4.2) 



(4.3) 



with v"" the components of the velocity field of the fluid. If the fluid is comoving with the expansion of the 
universe, i.e. , if it has no peculiar motion with respect to the general expansion, then w" = (1,0,0,0) in the 
coordinates (2.1) and the energy- momentum tensor is determined completely from the energy density, Tu — s 
and the trace T = —e + 3p. If the quantum state of the scalar field is spatially homogeneous and isotropic, then 
2($Vi$) == Vi(<I'^) = 0, and it foll ows immediately that the expectation value of its energy-momentum tensor is 
precisely of the perfect fluid form (43). Hence it suffices to consider 



{Tu)'^'e 



^(^'> - ^(*^3*> + 6Ci/($$) + 3C [h^ + ^) ($' 



($^ 



(r>'l'g'^''(rafc) = (6e-i) 



$2) _ ($A3$) + ei?($^ 



(6$-2)m2($2) 



(4.4) 



with 3p = 5'-' {Tij ) . Each of the expectation values of bilinears o f the scalar field in this formula are e asily computed 
as mode su ms (o r inte grals) by inserting the mode expansion (2.24), making use of the definition (2.3C), and the 
properties (^.21) and (^.23) of the harmonic functions. We find 
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(6^-1) 



2^2„3 / dA*(fc) gfclAf + ""^^^^^' I dfi{k)ak 



lAf 
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(4.5) 



with the isotropic pressure given by 



ie+{T)) 



(4.6) 



and 



d^(fc) 



J^k^dk K = 0, -1 



(4.7) 



The covariant conservation of the energy-momentum tensor may be checked by taking the time derivative of 
e and using the equations of motion (2.19) or (2.25). By exphcit calculation one can verify that 



e + 3H{e+p) = 0, 



(4.8) 



provided that in the form ( |4.5| ) one can take the time derivative inside the mode sums. Of course this calculation 
is still formal since the mode sums actually diverge and require the introduction of a cutoff or subtractions to 
become completely well-defined. In order to be ph ysic ally meaningful this regularization and renormalization 
procedure must preserve the conservation equation ( 4 



Now we may use the relations (3.1), (p.2|), and (3.7), to write the three quantities |/fc| , |/fe| , and Re(/^/fe) 



appearing in the energy density and trace in terms of the three bilinears Afk, Ti-k, and 2k and the two, as yet 
unspecified, functions Wk and Vk, which define the adiabatic basis. Explicitly, we have 



Cfc |/fc(i)| = — ;rrrrTT^ ^ 



2Wk{t) 



cTk \MtW 



2Wk{t) 



Wkit) 



Tlkit) 
Wk{t) 



-WliV) 



afeRe[/,*(t)/fc(t)]- 



'^^^'^V.(t) + ^y.(t)+x.(t). 



4W^fe(t) 



2Wk{t) 



V^{t) 



-Ik{t)Vk{t). 



(4.9) 



Hence we obtain 
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(4.10) 



where 
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(4.12) 
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2c.2 + 2iJ+4- 


2 y 


+ (6e- 


"1] 



+ (6C - 1) (^2^^,' - ^ + mVk 

Pf = ^(V^/. - i^) + ^(6e - l)m - Vk) . (4.13) 

We are now in a position to determine the specific choice of the functions Wk and Vk , by the requirement 
that all the power law divergences of the energy- momentum tensor should be contained in the vacuum zero point 
t erms , i.e., in the terms proportional to the factors of 1/2 multiplying e^, T^, and p^, respectively, in Eqn. 



(4.10). The adiabatic expansion of the mode functions fk may be employed to isolate these divergences. It has 



been shown by several authors p^,p5|,pq| that the fourth order adiabatic expansion of the frequency Wk{t)^ if used 



in (3.1) for the mode function and then substituted into the expressions for the energy- momentum tensor above, 
will reproduce all of the quartic, quadratic, and logarithmic divergences. The physical reason for this fact is that 
for any smoothl y va rying a{t) the large k behavior of the exact mode functions is determined by the WKB-like 
adiabatic form (^^) , and the local ultraviolet divergences of the energy- momentum tensor are identical to those 
obtained in the adiabatic expansion of the mode functions. Let us consider the matching of the large k behavior 
of the exact and adiabatic mode functions order by order in the adiabatic expansion. 

If we were to take Wk = Wf. = ojk and Vk = 0, then we would isolate the leading (i.e., quartic) divergences in 
the energy density and pressure. Indeed, to lowest adiabatic order we obtain the quartic and quadratic divergent 
terms. 



£ 



''' = ^j'^i^)-^ 



T'''--J^ f Mk)-. (4.14) 

47r^a'^ J ujk 

By differentiating e'"' with respect to comovin g ti me, one can easily show that this adiabatic order zero energy- 



momentum is formally conserved. It obeys (4.8) provided that the time derivative commutes with the sum 
(integral) over k, i.e., that any cutoff in the mode sum over comoving momentum is time independent. Thus 
subtracting it from the full expressions preserves the conservation law required by general coordinate invariance. 



even though the forms (4.14|) do not correspond to an explicitly covariant local counterterm in the effective action. 



The reason for this lack of manifest covariance is that the adiabatic expansion treats time differently than 
space, and in fact cutting off the sums in spatial momentum k corresponds precisely to a point separation in 
the spacelike hypersurface of con stant t [0-n5| . When this is recognized, then the correspondence between the 



adiabatic order zero expressions ( 4.14 ) and the subtraction of covariant countertems generated by point splitting 
in the spacelike direction may be made fully manifest and becomes completely justified. In the language of 
general renormalization theory, power law divergences are 'non-universal,' in the sense that their precise form 
depends on the regularization scheme employed, which in and of itself has no physical significance. This is 
clear from the fact that in dimensional regularization power law divergences do not appear at all. Hence any 
convenient scheme, covariant or not, may be employed to isolate and remove these divergences, although the 
form of the subtractions will not be explicitly covariant if the regularization scheme is not, and the existence of a 
covariant point splitting procedure equivalent to the non-covariant adiabatic subtraction is nec essary to establish 



the latter's validity. The fact that the subtraction of the non-covariant divergent sums such as (4.14) in adiabatic 



regularization nevertheless preserves the conservation equation (4.8) is a posteriori evidence enough that the 
procedure is completely consistent with general covariance, which is all that we require from a practical point of 
view. 



The subtraction of the adiabatic order zero expressions from the full energy density and pressure in (4.10) 



leaves subleading quadratic divergences still present. Hence we must match the function Wk and Vk to the 
adiabatic expansion to at least one higher order. Suppose that we choose to match Wk = Wj^ to second 
adiabatic order. We should allow then for Vk to enter at adiabatic order one (since VJ? and HVk, which are 
then adiabatic order two, appear in the energy-momentum tensor). Thus to this order we could choose Vk = 

— W^ /WJ^ = —iOk/uJk- Indeed this choice will match the quadratic divergences precisely, and one can verify 
that the second order expressions, 
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4^1 



2-^ 
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H^m" 



2H + iH^ + ^\- 






(4.15) 



Wi. /Wj. in order to match 



obey formal conservation on their own and may be subtracted from the fuU expressions in (4.1C). Again the fuU 
justification of this subtraction requires the covariant point sphtting analysis. At this order all the non- universal 
power law divergences are removed. 

Proceeding in this way one more time we might choose Wk — WJ: and Vk — —y}^ 
the remaining logarithmic divergences in the vacuum sector and leave behind completely finite state dependent 
terms in e and p. It is necessary to subtract up to adiabatic order four to obtain the correct finite trace anomaly 
in the conformally invariant limit. However, the remaining logarithmic divergences in the energy-momentum 
tensor (unlike the power divergences) correspond directly to local covariant counterterms in the quantum effective 
action, namely the fourth order local invariants, C abcd.C"'^'^''' and B? . The variation of the first of these vanishes 
in conformally flat FRW metrics, and a possible third linear combination of RabR""^ and B? vanishes by the 
invariance of the topological Gauss-Bonnet density, Rabcd.R"''"^'^ — ARabR"''' + R^, under local variations of the 
metric. He nce on ly the loga rithmic cutoff dependence of the coefficient of t he co variantly conserved tensor ^^^iJat, 
defined by (2.14) and (2.15), remains after the adiabatic order two terms (4.15) have been subtracted. In fact, it 
is easy to see that the fourth adiabatic energy density and trace are given by 



sb/^"»'i!(«-i>'"*' 



{H^ - 2HH - 6HH^ - H"" 




W-Di'^^ 



6HH^ + 2HH -H^ - —H' 



^-ii/^-'^if («-['""' 
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tH 
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(4.16) 



^ ( ^ - 7 ) ( -J + 4-T^ + 2-^i?^ + lei/" + 30HH^ + AHH + 3H^ 



As is evident from these expressions the adiabatic terms of order four are of two kinds: logarithmically divergent 
pieces, which arc proportional to the local geometric tensor ^^' Hat, and finite pieces which are given in terms of the 
quantities a, H, H, H, and H- The logarithmic cutoff dependence of (Tab) may be absorbed into a renormalization 
of the coefficient of the tensor ^^'Hab in the backreaction equations, i.e., by a renormalization of the coupling 
constant of the local R^ term in the effective action ||9[. Since such a term has to be introduced in principle 
into the backreaction equations in any case, we may simply include an explicit logarithmic cutoff dependence 
(with the correct coefficient) in the bare R? coupling to cancel the remaining logarithmic cutoff dependence in 
the vacuum energy-momentum tensor of the quantum field. Th e re sulting backreaction equations will be fully 
cutoff independent and obey the covariant conservation relation (4.8). 

In electrodynamic backreaction problems the analogous procedure has been employed in practical numerical 
implementations p9| . The logarithmic divergence in the electric current expectation value (j") in that case 
does not need to be removed by explicit subtractions. Rather the logarithmic cutoff dependence of the bare 
electric charge, cq, is used to cancel the cutoff dependence of the expectation value of the current {]"') in the 
full backreaction equations, {1 / e^)dbF°-'^ = (j°). Provided that the momentum cutoff is chosen much larger than 
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the inverse scale of temporal variation of the electric current and fields F"^'', it can be demonstrated numerically 
that the results are independent of the cutoff, as they should be. This can be checked explicitly by running the 
numerical code with several different cutoffs and verifying that the evolution is unchanged, provided that the 
bare charge is rescaled logarithmically with the changed cutoff in such a way as to keep the renormalized charge 
fixed. In other words, rather than striving to eliminate it explicitly, the remaining logarithmic cutoff dependence 
in the expectation value of Tab may actually be used to our advantage, in order to check the cutoff independence 
in the numerical evolution of the backreaction equations. 

This procedure may be implemented in the semiclassical backreaction equations as follows. We propose to 
avoid explicit subtraction of the logarithmic cutoff dependence of the energy and pressure, by checking that it 
cancels against the logarithmic cutoff dependence of the bare fourth order coupling multiplying the tensor ^^'Hab, 
which appears in the semiclassical Einstein equation of backreaction. If we denote the bare coupling of the i?^ 
term in the action by l/ofQ then the logarithmic cutoff dependence in the backreaction equations will be removed 
if 

where A is the cutoff in the comoving momentum sum (integral) and /i is an arbitrary finite renormalization 
scale. The finite adiabatic order four pieces in the energy-momentum tensor are not taken into account by this 
procedure, so they can be added back in by hand to the backreaction equations. 

Besides being conceptually and technically somewhat simpler than the usual approach of explicitly subtracting 
all terms up to adiabatic order four, this proposal is a significant simplification for practical implementation of the 
backreaction equations on a computer since WJ^ as given by (3.6) involves both three and four time derivatives of 



the scale factor a{t) |30| . These higher order derivatives present problems for the standard approach to numerical 
backreaction calculations, when they occur inside the mode sums, since one cannot determine the form of (Tab) to 
step the equations forward in time without first knowing these higher derivatives, which are themselves determined 
through the backreaction equations by (Tab)- 

Since in the general FRW backreaction problem we will not remove the logarithmic cutoff dependence in the 
energy and pressure connected with the geometric tensor ^^'Hab, we will not need to match the functions Wk 
and Vfc to their full fourth adiabatic order values. However even if these functions are matched only to second 
adiabatic order, fourth order terms in the 'vacuum' contributions to the energy and pressure, e^ and p^, will be 
generated inevitably, from the square of adiabatic order two terms in Wk- In order to cancel these terms in e^ 

and p^ we will choose Wk to be W^ plus those additional terms on the first line of (3J) which fall off only as 



f ast as k~'^ at large k. Likewise we will choose Vk to be —ujk/oJk plus the one additional term on the first line of 
( |3.5| ) which falls off only as fast as k~^ at large k. 

Since we are not attempting to isolate the adiabatic order four logarithmic divergences, and wish only to 

(2) 

avoid introducing spurious fourth order divergences from the square of Wf^ , we may evaluate the fourth order 
terms in Wk and Vk on a spacetime with vanishing ^^'Hab- This is the minimal choice of Wk and Vk, which 
leaves behind only a logarithmic cutoff dependence in e and p proportional to the geometric tensor ^^' Hab- The 



vanishing of ^^^ Hab, with components given in eq. (2.15), allows us to make the following replacements 
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R -^ -3RH -^ --R IH - —] , (4.18) 

2 V a^/ 

when evaluating Wk and Vk- Hence the higher derivative terms in the first line of (|3.6[) become 



and we finally define: 



R + 5RH + 2RH + 6RH^ ^ ^ [R-^] , (4.19) 
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(4.20) 



(4.21) 



to specify fully the adiabatic particle number basis for the scalar field of arbitrary mass and curvature coupling 
in a general FRW spacetime. 

Some comments about this definition and the corresponding choice of basis bear emphasizing. First, by 
construction these choices for Wk and Vk match all the quartic and quadratic power law divergences in the energy 
and pressure in a general FRW spacetime. Thus the subtraction of the adiabatic two expressions (4.15) from the 
bare energy density and trace will leave behind only the logarithmic cutoff dependence proportional to ^^'Hat, 
whenever it is non-vanishing. 

Second, with these requirements, the choice of Wk and Vk, and hence the adiabatic basis they define is the 
almost unique, minimal choice. Since the adiabatic (order zero and two) matching is required to eliminate the 
non-universal power divergences in the energy-momentum tensor, and since the matching of any higher adiabatic 
orders requires higher order time derivatives of the scale factor which we exclude, the only freedom left in the 
choice of Wk and Vk are adiabatic order four terms in Wk and third order terms in Vk, which either give rise 
to local '^'^'^Hah terms in {Tah) or fall off faster than 1/k^ and hence lead to finite reapportionment of the terms 
in the energy-momentum tensor between 'particle' and 'vacuum' contributions. This residual ambiguity in the 
separation of the energy-momentum terms into particle and vacuum contributions cannot be removed except by 
an essentially arbitrary choice. It reflects the necessary uncertainty in the definition of the 'particle' concept 
in time varying external fields. However the residual ambiguity is fairly mild and does not affect severely the 
physical interpretation of particles in semiclassical cos molo gy where the variations of the metric are assumed 
small in Planck units. Our definitions of Wk and Vk in (4.20) and (4.21) are the minimal ones, involving no more 
than second derivatives of the metric, that make such an interpretation of the separation of {Tab) into quantum 
vacuum and quasiclassical particle components possible for slowly varying a{t). 

Third, for the massless, conformally coupled scalar, m — and ^ — ^, th e definitions ( 4.2C ) and ( 4.21 ) 
simplify considerably: Wk = k/a and Vk = H. Referring back to (2.26) and ( |3.18 ) we observe that the adiabatic 
particle number A4 is constant in this case, i.e., A/fe = Nk. Hence there is no creation of massless, conformally 
coupled field quanta in arbitrary FRW spacetimes, as expected from conformal invariance. 

Fourth, for the massless, minimally coupled scalar field, ?7^ = ^ = 0, 
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(4.22) 



We observe that Wk is everywhere non-negative for R > 0, vanishing only for the spatially constant mode in 
the case of closed spatial sections, k = k = 1, and then only at times for which R — 12 /a^. With these special 
exclusions the definition of adiabatic particle basis specified by (4.20) and (4.21) is meaningful even in this quite 
infrared sensitive case, which is the most analogous to linearized graviton fiuctuations. 

In this paper we restrict ourselves to a test field approximation, and do not consider the backreaction of the 
created 'particles' in the background geometry. For the special case of de Sitter, which will be the particular 
spacetime under study, the tensor ^^'Hab vanishes. Hence in this fixed background the one remaining logarith- 
mic cutoff dependence proportional to the geometric tensor ^^' Hab (after subtraction of the adiabatic order two 
expressions) will vanish identically, and the resulting renormalized energy-momentum tensor is fully cutoff inde- 
pendent, without any additional countertems. In order to compare the value of the finite (Tab) with previous 
authors, we perform the final adiabatic subtraction of order four as well. The fact that ^^•'ff qfc is zero fo r de Sitter 
simplifies this order four subtraction, so that we only need to subtract the finite pieces of ( 4.16| ) and (4.16), and 
the logarithmic cutoff dependent counterterm is not needed to obtain the renormalized energy-momentum tensor. 

Examining in detail the adiabatic 'va cuum ' component of the energy and pressure, obtained by subtraction 
of the second order adiabatic expressions (|4.15| ) , 
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(4.25) 



and using the definitions of Wk and Vk in ( 4.2C ) and ( 4.21 ), we ob serve that the sums (integrals) for the these 
vacuum contributions to the energy density and pressure in (4.23) are exphcitly convergent for arbitrary mass 
and curvature coupHng ^ in a general FRW spacetime. Hence the remaining logarith mic cutoff dependence of 
{Tab) proportional to ^^' Hat must reside only in the terms linear in N'k, "R-k, and Ik in (4.10), and this remaining 
logarithmic cutoff dependence vanishes as well in the special case of de Sitter spacetime. 

An additional technical point is the fact that when performing the adiabatic subtraction one should use the 
continuous measure, and not the discrete one, even in the case of a closed cosmological model [Q. This scheme 
is necessary to compute finite Casimir energies and the correct trace anomaly, for example. We only point out 
at this stage the relevant steps needed to carry out this subtraction. The details of this calculation can be found 
in|I|. 

The terms needed are called the Plana terms, and are given by 
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(4.27) 



We need to calculate the Plana terms for n = — 1, n = 1, and n = 3, as these values of n correspond to the 
divergences present in the theory. It is easy to see that the calculation of the Plana term for n = 3 requires the 
introduction of an infrared regulator (. With this is mind, we have 



Vi 



1 


V 1 = 


-1 


12U 




12 



V-3 = log C + 7£ = logC - ^(1) = log C + 0.5772157 . 



(4.28) 



The infrared regulator drops out of the final answer for the renormalized (Tab) when the state is IR finite. We 
present results for the renormalized (Tab) for the special case of a de Sitter universe in the next section, with 
these additional finite Plana terms subtracted (as described in |Q) in order to compare our results to the earher 
literature. 
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V. PARTICLE CREATION IN DE SITTER SPACETIME 



The line element of de Sitter space with k = +1 is given by (2J) and (2^) with 

O'desi't) — H"^ cosh Ht . 



(5.1) 



Because of its higher degree of symmetry and special status both in the preceeding discussion of renormalization 
of (Tab) and in cosmological models of the early universe, we will apply our decomposition of the renormalized 
energy-momentum tensor of the scalar field first to de Sitter spacetime. The de Sitter invariant 'vacuum' state (or 
variants thereof) is by far the most commonly discussed in the literature pi| |. However, from the point of view of 
the initial value problem for the scalar field theory, specified on a complete spacelike Cauchy surface, there is no 
a priori reason to impose the global 50(4, 1) invariance of de Sitter spacetime. Our general framework permits 
us to consider any spatially homogeneous and isotropic initial state consistent with finite, renormalized energy 
and pressure. 

This general initial state of the scalar field i s spe cified by initial data on the mode functions, /fc(io) smd fkitg), 
constrained to obey the Wronskian condition ( 2.29 ), together with the initial particle number density in phase 
space Afkito)- By convention and without any loss of generality we will retain our previous definition of Nk t o 
be the initial particle number density with respect to the initial adiabatic mode basis specified by (3J_), ( 4.20 ), 
and ( 4.2l[ ). This means that for general initial data we must allow /fc(to) and fkita) to differ from their adiabatic 

vacuum valu es fk{to) and fk{to), respectively. Consequently, the Bogohubov coefficients ak{to) and (ik{to) do not 
satisfy ( ^.14D and Mk{to) ^ Nk in general. 

The global phase freedom in the mode functions allows us to choose /fc(io) real. Thus we may define the two 
real functions of fc, Wk and Vk, by 



fkito) = 
fkito) = 



'iwk + y ) fkito), 



(5.2) 



which is the general form of initial data satisfying the Wronkskian condition (2.29). The physical condition that 
the initial state have finite energy and pressure requires that Wk and Vk must match the adiabatic vacuum values 
Wk and Vk, defined by (4.20) and (4.21), up to terms that fall off sufficiently fast at large k, i.e.. 



\wk-Wk\ <0{k-'') 
Vk\<0(k-^) 



\Vk 



(5.3) 



as A: — > cx). Likewise we must require that the initial particle distribution have finite energy and pressure, so that 



Nk<0 (k 



-3\ 



as k 



(5.4) 



Any choice of finite initial data in the form of the three functions Wk, Vk, and Nk satisfying these conditions at 
large k have finite initial renormalized energy and pressure and are physically allowed. This condition of finite 
initial energy and pressure is sufficient to guarantee finite energy and pressure for the scalar field at all subsequent 
times. 

In the special case of spacetimes for which ^^' Hat — 0, such as de Sitter space, we may saturate the inequalities 
( ^.3[ ) and (5.4), since the logarithmic ultraviolet divergences in (Tab) this would lead to in a general FRW space 
are proportional to ^^' Hab and vanish when ^^' Hab vanishes. 

The numerical solution of the initial value problem proceeds as follows. We first define initial conditions as 
given in eq. (5.2), and then solve the mode equations using a sixth order Runge-Kutta integrator. All quantities 
of interest can be derived from the numerically evaluated fk(t) and fk{t)- Computations of the energy- momentum 
tensor involve direct summations (not integrals) since we are considering closed spatial sections. Since all the 
modes are independent the calculation is well-suited to a parallel computer. 

The numerical scheme preserves the covariant conservation of both the bare and the renormalized {Tab)- 
The conservation of the bare tensor is a direct result of the equation of motion for the field $(i, x). Because the 
adiabatic subtraction procedure at any given order preserves the covariant conservation of the renormalized {Tab), 
{Tab) R must be conserved as well. In Fig. |l| we show explicitly the covariant conservation of the renormalized 
energy-momentum tensor for to = 1.0, ^ = 1/6, t^ = 2.0, and a momentum cutoff of A = 400. 



16 



1 1 1 r 



dEnergy/dt 




3 H(Energy+ Pressure) 



3.6 



3.8 



FIG. 1. The top line plots ^^, the bottom line plots 3//(eR +pr), and the sum of both is vanishing, showing that the 
renormalized energy-momentum tensor is covariantly conserved. 

One can also show explicitly that our results are cutoff independent. That is, if we run with a momentum 
cutoff A (the total number of modes is A) , all renormalized quantities remain cutoff independent until we reach 
the time at which all the modes have crossed the horizon. That is, we can trust our results until times of order 
of iA, with ^A given by cosh(ffiA) ~ A. 

In Figs. 2 and 3 we plot the particle number Afk{t) in de Sitter space for various values of fc, m, and ^, 
starting at to = ~8 and to — (in units oi H ~ -R/12) with adiabatic vacuum initial conditions, and the choices 
Wk = Wk and Vk = Vk- When m and ^ are such that 



7 



M^ 
W 



9 

4 



ni 



12^ 



>0, 



(5.5) 



the particle number increases sharply (with a rise time of order H ^) when the physical wavelength of the mode 
crosses the de Sitter horizon, i.e., when t = tr, where 



cos\v{Htc) — k . 



(5.6) 



When the initial time to ^ 0, the de Sitter universe (in closed spatial coordinates) is always expanding thereafter 
and each mode can cross the horizon and go through its particle creation event at most once. This situation is 
shown in Fig. || and Fig. ^ If the initial time to < 0, the spatial sections first contract, then expand, and some 
modes can c ross the horizon twice, once on the way in and then again on the way out, corresponding to the two 
solutions of (5.6) at t = ±tc. These modes then go through two separate bursts o f pa rticle creation, as illustrated 
in Fig. 0. The behavior of the horizon crossing time with k in accordance with (|5.6|) is shown in Fig. ra and Fig. 
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FIG. 2. Particle number for m — 0.6, ^ — 1/6, to = ~8, and different values of comoving momontuni, from left to right 
k = 1,2,10,20,25,50,100. 
0.005 




FIG. 3. Particle number for m=1.0, ^=1/6, to = 0, and different values of comoving momentum, from left to right k—1, 
20, 50, 100, 200, 250, 500, 1000. 
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10 

FIG. 4. Particle number for m=1.6, ^=0, to = 0, and different values of comoving momentum, from left to right, k—1, 
20, 50, 100, 200, 250, 500, 1000. 

In Figs. 5 and 6 we plot the logarithms of the comoving momentum as a function of their crossing time i;.- 
The lines correspond to the fits, O.Sexpt and 0.16 expi, for ^ = 1/6 and ^ = 0, respectively, in agreement with 
the previous considerations regarding the nature of the particle production at horizon crossing and earher work 
on stochastic inflation [p2|,^. 



o 



1000 



100 



10 



1 







1 1 


1 


1 




1 




































- 












^-' 


- 


















- 












,o' 


- 


- 










^^ 


''' 


- 












J^" 














__, 


^iS> 














,^' 








































































- 






^Sf' 








- 


■ 






•' 








' 


- 




^^ 










- 






.'*'' 












































- 




^^ 










- 


































- 


^-"^ 












- 




,.■'' 
















^,''' 
















o 


















1 1 


1 


1 




1 





1 23456789 

t_c 

FIG. 5. Momentum versus comoving crossing time for m=1.0, 5=1/6, and fc=l, 20, 50, 100, 200, 250, 500, 1000. 
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FIG. 6. Momentum versus comoving crossing time for m = 1.6, ^ = 0, and fc=l, 20, 50, 100, 200, 250, 500, 1000. 

This corresponds to the physical picture of particle creation as a parametric resonance between the wave 
function of the particle mode and the geometry, which is maximized when the physical wavelength of the mode 
is the same as the horizon scale. 

In Figs. 7 and 8 we have plotted the renormalized values of the energy, trace, and pressure for the cases 
considered above. 
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FIG. 7. Energy density, trace, and pressure for to = 0, m, = 1.0, and ^ = 1/6. 
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FIG. 8. Energy density, trace, and pressure for to = 0, m=1.6, and ^=0. 

The most important feature of the figures to note is that the particle number goes to a constant which is 
almost independent of k at late times, when the modes freezes out. The value of the constant plateau value of 
Afe at late times and large k can be compared with that calculated by one of us p3 analytically some years 
ago, namely Afk — > sinh~ (TT-f), assuming an initial 'in' vacuum state corresponding to pure positive frequency 
as to — > — oo. Although this 'in' state is a higher order adiabatic vacuum state than the one considered in the 
present work, the difference between the two goes to zero rapidly for large 7. We remark that if we had chosen 
a different separation of particle and vacuum components of (Tab) by making a slightly different definition of Wk 
and Vk at adiabatic order four, then the detailed time profile of the Afk near the horizon crossing would be slightly 
different, but the location of the crossing time and the plateau values of particle number at later times observed 
in the figures would not be altered. 

The fact that the particle number density in any given k mode does not redshift away despite the exponential 
expansion is related to the fact that these created particles do not satisfy a standard dust or matter equation of 
state with e > 0, p > 0. Rather as Figs, uj and show, the matter plus vacuum contributions to the renormalized 
energy and pressure satisfy the de Sitter equation of state, e — —p at late times. In this free field theory there is 
no scattering between the created particles, and a 'normal' equation of state cannot be established. The form of 
the equation of state is expected to change as soon as interactions are introduced, no matter how weak, since the 
particles have a finite density and an unlimited interval of time to interact. 

Next we consider the massless cases {m — 0) for both conformal (^ — 1/6) and minimal (^ = 0) couplings. 
As already remarked, in the conformal case there is no particle creation at all, so the energy and pressure remain 
strictly constant and agree with that of the Bunch-Davies vacuum pi[ given by 



{Tal 



\B-D 



o,e 



gab R 
64^2160 



i/4 



9607r- 



;gab 



(5.7) 



which corresponds to the value of the anomalous trace. 

The results for the massless, minimally coupled field are shown in Figs. || and |l^. We point out the fact 
that the renormalized value of {^'^)r grows linearly in time with a slope of -^, that is independent of the initial 
conditions ||23,|3l. 
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FIG. 9. Renormalized value of {$^)b as a function of comoving time for m = 0, ^ — 0, and to = 1.0. 
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FIG. 10. Energy, trace, and pressure as functions of comoving time for m = 0, ^ = 0, and to = 1-0. 

Notice that for the massless, minimally coupled case the late time behavior of the renormalized energy, trace, 
and pressure are determined by the value given in reference p2[ , and it does not go to the value calculated in 
reference [£1[. This difference can be understood as a finite additional term that must be added to the Bunch- 
Davies result when the spatially homogeneous k = 1 mode is handled properly. 

It is interesting to study the case of the exactly massless, minimally coupled scalar field for a wider variety 
of initial conditions. If a non-zero number of particles Nk in a finite number of modes is added, the asymptotic 
values of the renormalized (Tab) still approaches the de Sitter invariant value [p2[ 



{Tab) a 



nar\ Tdab : 



(5.8) 



even though the initial state and transient behavior in those cases (non-vacuum initial conditions) are not de 
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Sitter invariant. 

The fact that in all the above cases, even for de Sitter non-invariant initial conditions and even for the 
massless, minimally coupled field, the total renormalized energy-momentum tensor goes to the de Sitter form, 
£ — > —p — > constant, is an important general conclusion of our study. Qualitatively, one can understand this 
result in the following way. At late times, t 3> H~^, a{t) grows exponentially in de Sitter spacetime, and a 
very large number of modes have wavenumber k ^ a, i.e., physical wavelengths that have been redshifted out 
of the de Sitter horizon and frozen out. Those wavenumbers for which k is still greater than a(t) behave like 
adiabatic vacuum modes and are removed by the subtraction procedure. Hence the sums (integrals) over k in the 
renormalized energy- momentum tensor are cut off at fcmax of order a{t). Since the initial conditions differ from the 
de Sitter invariant vacuum data only for some finite set of wavenumbers k < ko (say), these small wavenumbers 
make a decreasingly small contribution to the mode sum, which is controlled by the wide band of wavenumbers 
for which feg ^ fc < a{t). The summand (integrand) in this wide k band appearing in the energy density, pressure 
or trace is a function only of physical wavenumber k/a (since it is independent of the initial data). Examining 
the behavior of the mode functions for late times and expanding the integrand in e or p in powers of physical 
momentum it is easy to see that 

Hence all terms in the summand which do not grow faster than n — 3 for large a{t) yield a constant contribution 
to (Tab) in the late time limit, independent of the initial data. The largest power n which can occur in the 
mode sums for s and p is controlled by the behavior of the mode functions fk at late times. Referring back to 
Eqs. (2.25|) and (|2.26D it is clear that the late time behavior of fk{t) is oscillatory, exp{±i^Ht) for 7^ > but 



exponential, exp(±|7|iJi) for 7^ < 0. Since (Tat) is bilinear in the mode functions the power n can be as large 
as 2I7I for 7 pure imaginary. Provided 7^ > — |, all n < 3, and the exponential behavior cannot overcome the 
normal kinematic redshift a^^, there are no extreme infrared singular terms in the summand, and we conclude 



from (5.9) that both s and p should tend asymptotically to a constant, independent of the initial data at late 
times. By the covariant conservation equation ( |4.§| ) for the renormalized (Tab) this implies the de Sitter equation 
of state e -^ —p at late times, which is precisely what the numerical results show. 



Conversely, if 7^ < — | , we should expect the integrals for the finite renormalized energy density and pressure 
of the scalar field to contain terms with n = 2\j\ > 3 and to grow with the scale factor like a^''^''^^ [Q. Details of 
this case and the analytic calculation of this behavior will be presented elsewhere [M . 

The case 7^ = — | is marginal. For this value of 7^ the power n = 3 is possible (but not required) in the mode 
sum. If present, it leads to a logarithm, log a — > Ht and linear growth in cosmic time for the components of the 
energy-momentum tensor. One has to examine the energy-momentum tensor in more detail to determine if such 



logarithmic behavior exists or not, as a general argument based on the asymptotics of the mode equation (2.25) 
will not suffice. In the adiabatic 'vacuum' and matter contributions separately we observe this linear growth for 
the massless, minimally coupled scalar. This beh avior of the 'vacuum' component is easily understood from the 
asymptotic forms of the functions Wk and Vk in ( 4.22| ) in the massless, minimally coupled case, since they imply 
that 



1 1, 



iff . (5.0, 



for k ^ a. This is the marginal n — 3 power behavior which leads to a linear growth (for the energy density 
of the vacuum contribution) in comoving time with slope (Stt^)"^. However, it turns out that for the exactly 
massless, minimally coupled case {m^ = ^ = 0) the complete, vacuum plus particle, renormalized energy and 
pressure do not have such n = 3 logarithmic divergences, as our numerical studies show. In this case there is an 
exactly compensating growth in the matter contribution which cancels the vacuum piece. This cancellation is a 
consequence of the exact field translational symmetry (<& ^ $ + constant ) in the massless, minimally coupled 
case which prevents a term proportional to ($^) from appearing in (Tab). 

Next we consider 7^ = — |, i.e., m^ + 12^H^ = 0, with m^ and ^ separately different from zero, which has 
t he same mode equation but no such symmetry under field translations. Because the energy-momentum tensor 



(4.4) depends on m^ and ^ in distinctly different ways, in these cases there is a generic term proportional to ($^) 
in s and p, and hence the exact cancellation of the linear growth in t between vacuum and particle contributions 
no longer occurs. The total energy-momentum tensor now grows logarithmically with a{t) or linearly with t in de 
Sitter space. This behavior is illustrated in Fig. |Tl|. 
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FIG. 11. Energy, trace, and pressure as functions of comoving time for m = 0.3, ^ = —0.0075, and to = 1.0. 
The coefficients of tliis linear growth can be determined analytically from the known behavior of ($^), namely 



(T>- 



47r2 
' 47r2_ 

47r2 



^ + £0, 
t + To, 



with the constants Eq and po constrained by the conservation equation ( |4.S| ) to obey 



and the constant Tq determined by the equation 



Tq = 3po - £o ■ 



(5.11) 



(5.12) 



(5.13) 



Hence for ^ > (and m^ < 0) this linear growth decreases the effective cosmological constant over time, whereas 
for 'C < (and m^ > 0) it increases it. In either case the backreaction of the scalar quantum field cannot be 
neglected at late times. It remains to be seen whether this linearly growing behavior carries over to the physically 
more relevant case of one-loop metric fluctuations themselves, and if so with which sign. 

Finally, it is interesting to study the case of minimal coupling and non- vanishing but small real mass. In this 
case there are two time scales: a short time scale on which the AUen-Folacci vacuum is approached, and a longer 
time scale controlled by 1/m on which there is a smooth transition to the Bunch-Davies vacuum. Details of this 
case and a reliable analytic calculation of the behavior observed will be presented elsewhere 124] . 

VI. CONCLUSIONS 

We have presented a general method to study the initial value problem of a quantum field in a FRW spacetime. 
The two basic ingredients are: an initial adiabatic state (of order four), and an adiabatic number basis, with respect 
to which we define a particle number. This basis allows one to write the bare energy-momentum tensor as a sum 
of two contributions: a vacuum polarization term, and a particle production term. The advantage of the basis 
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presented here is that the non-universal power law divergences of {Tab) are isolated in the vacuum polarization 
piece, and the matter contribution displays physically reasonable behavior, as evidenced by the particle creation 
process in de Sitter space. We have implemented a numerical scheme incorporating this theoretical approach 
which preserves the covariant conservation of (Tab), and yields a renormalized (Tab) which is cutoff independent. 

The numerical results for a de Sitter background can be summarized as follows. For the case of large mass 
(7 ^ 0) particle production takes place at horizon crossing, and the late time behavior of the renormalized energy- 
momentum tensor agrees with that of the Bunch-Davies vacuum. In the conformally invariant case {m = 0, ^ = 
1/6) there is no particle production, and the renormalized value of the energy-momentum tensor corresponds to 
the well-known trace anomaly. The massless minimally coupled field (m = 0, ^ = 0) appears to be a marginal case. 
Non-de Sitter invariant initial conditions are driven to the AUen-Folacci vacuum, and not to the Bunch-Davies 
vacuum. In none of the above cases does (Tab) vanish for late times, that is, it is not red-shifted away due to the 
expansion of the universe. Backreaction effects, though not required by this behavior cannot be ruled out when 
higher order corrections due to self-interactions or fluctuations in the background metric are included. Finally, for 
m^ + ^R = 0, but m^ and ^ individually different from zero, the total energy density, trace, and pressure of a free 
scalar field grow without bound, logarithmically with the scale factor. In this case it is certain that backreaction 
effects must be included self-consistently, since the energy density and pressure must eventually grow large enough 
to influence the background geometry and invalidate the test field approximation. 

We believe that these results merit a fuller investigation of the effects of self-interactions and quantized 
gravitational fluctuations in a complete backreaction calculation. With respect to the latter, an interesting 
extension of the present study of free scalar field theory would be to the second order energy-momentum tensor 
of quantized metric fiuctuations themselves. Barring a special cancellation which occurs in the ^ = ra = scalar 
case, we may expect that this energy-momentum tensor would contain terms which grow logarithmically with 
the scale factor, and dynamically infiuence the evolution of the mean geometry away from de Sitter space, as 
has been suggested by general arguments some time ago |35|] . The generalization to fully self- interacting fields 
is another interesting direction to pursue, since one would again expect the interactions to alter the de Sitter 
equation of state found in our study of the asymptotic behavior of massive fields. The formalism presented in 
this paper is intended to serve as a foundation from which the detailed analyses of these issues can begin, and 
the general backreaction problem in cosmological spacetimes with arbitrary finite adiabatic initial data can be 
posed and solved numerically. These extensions, applications to specific cosmological models, and the generation 
of density perturbations in the early universe are under current study. 
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APPENDIX A: ADIABATIC EXPANSION IN CONFORMAL AND COMOVING TIME 

Because of the loss of explicit covariance under time reparameterizations in the adiabatic regularization 
procedure, we establish in this Appendix the relationship between the adiabatic expansions in two different, 
commonly used time coordinates, namely conformal time 77 and the comoving time t employed in the present 
work. We will show that the two adiabatic expansions (comoving and conformal) give identical results for physical 
quantities such as (Tab) to any order in the expansion. 

The relation between the two time coordinates is given by ( |2.3|) of the text. We denote a derivative with 
respect to comoving time by a dot, and a derivative with respect to conformal time by a prime. If we define the 
conformal mode functions by the relation, 

9k{v)^a--2{t)fkit), (Al) 



then the mode equation ( 2.25 ) in comoving time implies the mode equation, 

9kiv) + ^Uv)9k{v)^0, (A2) 

with 
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mlirj) = C{v)u;l +{C-l) C{^)R{^) = a^^i + a^ (^ + ^ 



(A3) 



in confornial time. Hence the zeroth order adiabatic frequency for the adiabatic expansion in confornial time is 
a times the zeroth adiabatic frequency for the adiabatic expansion in comoving time, LUk- We wiU now show that 
this feature holds to any adiabatic order. 



Since the mode equations, (A2) or (2.25), both have the standard form of a time dependent harmonic osciUator 
equation, the adiabatic ansatz. 



9kiv) 



[2T,(r?)]i/2 
leads to an equation for Tk{ri) in confornial time, 



exp 



i dri'Tkiv') 



(A4) 



^2 „2 
1 I. = 117 



3 1 /dT 



1 1 d^Tfc 

which is exactly of the same form as that in comoving time, upon replacing Tfc by Wk, and tuk by flk, i.e., 

^'-^^''^AW^y dt ) 2Wu di2 ■ 



To zeroth adiabatic order we have established already that 

Let us now assume that this result is true at adiabatic order N — 1, i.e., 

and iterate to order N. Since 



dT 



(N-l) 



dr] 



d^Tf-^) 3 



iJW^f-i)+H.f-i) 



2i/2^f -1) + Hwt''-'^ + 3iJM^f ^^) + VFf -^) 



d?7^ 

substitution into the right side of ([A5|) and using ([A3|) gives the desired iteration to order N, namely, 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



,(7V)2 



o 3 1 



dT 



(N-l)' 



1 1 d^T 



2^iN-l) 



4 -^(^^-1)2 I dr/ 



2 T^^"^' ^V^ 



nl + ^ 



1 



dW^ 



(N-l)' 



1 1 d^W, 



2w(W-l) 



" 4 u/(^-i)2 



w, 



dt 



2w(A^-i) 



di2 



= a'wi''^' 



(AlO) 



Hence by induction we have proven that to any adiabatic order 

TW(,)=a(t)W^f)(i). 



(All) 



Moreover it is clear that this argument has relied only on the form of the mode equation (2.25) or ( A2), and their 
corresponding WKB definition of the frequencies Wk and T^. Since under an arbitrary differentiable transforma- 
tion of time coordinates, 



dt = b{t')dt' 



(A12) 
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the mode functions can be rescaled, 

fk{t)^b-Ht')gk{t'), (A13) 

in order for the first t' derivative terms in fk to be ehminated, the mode equation in the new time variable can 
be put in the standard form, 

gUt')+r^l{t')9,{t')^0, with ^2(^,) ^ ^2^2 _ 3 /Q ^ b^ ^ (^14) 



where now the prime denotes differentiation with respect to t' . Once in this form the adiabatic ansatz (A_5) may 
be inserted and the previous inductive argument may be repeated to yield 

T[^\t') = bWJi'^Ht) . (A15) 

for the arbitrary reparameterization of the time coordinate, t -^ t'{t). 

An immediate corollary of this result is that any quantity such as the energy density and pressure, developed 
up to a given adiabatic order N in the comoving time coordinate t corresponds exactly to the same quantity 
developed up to the same adiabatic order N in any other time coordinate t' . Since Tab transforms as a tensor, 
the energy density and pressure transform as 

£ ^ e' = &^e , 

p ^ p' = p , (AI6) 



under (A12). Hence up to the overall factor of b^ in the first of these relations, the energy density and pressure 
in either time coordinate are numerically equal at any order of the adiabatic expansion, and the renormalization 
procedure of subtracting the expansion up to adiabatic order two or four from the bare quantities is the same 
procedure in any time coordinate. The case t' = rj and b = a, relating conformal time to comoving time, is just 
a special case of this much more general result, so that the conformal adiabatic expansion which has been given, 
f or e xamp le, i n |ljJ2q] is exactly equivalent to the corresponding expansion in comoving time given by Eqns. 



(3.4) and (3.6) of the text, as may be checked by explicit comparison of the two sets of expressions. Hence the 



subtraction of the adiabatic order two expressions (4.15) from the bare energy density and pressure corresponds 



exactly to the same subtraction in conformal (or any other) time coordinate which leaves the spatial sections 
unchanged. Finally, since the fourth order logarithmic cutoff dependence of {Tab) is proportional to the geometric 
tensor ^^'Hab, this remaining cutoff dependence and the argument (given in the text) for how to handle it in 
backreaction calculations, applies equally well in any time coordinate. 
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